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EXPERIMENTAL INVESTIGATION OF FREE-CONVECTION HEAT TRANSFER IN VERTICAL TUBE 

AT LARGE GRASHOF NUMBERS ^ 

By K. H. Ci. Eckert and A. J. Diaguila 


SUMMARY 

An investigation was conducted to study free-con vection heat 
transfer in a stationary vertical tube closed at the bottom. The 
walls of the tube were heated^ and heated air in the tube was con- 
tinuously replaced by fresh cool air at the top. The tube was 
designed to provide a gravitational field with Grashof numbers 
of a magnitude comparable with those generated by the centrifu- 
gal field in rotating-blade coolant passages {10^ to 10^^). Local 
heat-transfer coefficients in the turbulent-flow range and the 
temperature field within the fluid were obtained. 

I leaf -transfer coefficients, converted to dimensionless Nusselt 
numbers, were correlated as a function of Grashof and Prandtl 
number parameters. When the heat-transfer coefficients were 
based on the temperature diflference in the horizontal plane 
under consideration, the Nusselt numbers in the turbulent range 
investigated averaged 35 percent below the known relations for 
vertical flat plates. 

Several tube-cover configurations that provi/led different 
methods of removing and reffiacing the heated air were investi- 
gated. The greatest drop in the temperature of the air in the 
tube occurred very near the tube wall. This drop indicated that 
a thin boundary layer existed along the heated wall surface. 
The temperature of the air in the central axis of the tube in- 
creased about 30^ F with distance from the air inlet. When 
the bottom of the tube was closed, the location of warmer lighter 
air below cooler heavier air created unstable conditions in the 
flow. The heat-transfer coefficients determined by use of local 
temperature differences between the wall and the air in the tube 
axis as well as the magnitude of the temperature increase along 
the tube axis were not influenced by the configuration changes. 

INTRODUCTION 

Interest in the study of free-conveetion heat transfer in- 
creased when it was realized that the centrifugal field in ro- 
tating components of propulsion systems often creates very 
strong free-convection currents. The existence of such 
Hows may, for example, influence considerably the heat 
transfer in air-cooled rotor blades of gas turbines and account 
for certain observations in the combustion chambers of ro- 
tating ram jets used to propel the rotor blades of helicojitors. 
Furthermore, it was found that in liquid-cooled turbine rotor 
blades the fi*ee-convection currents alone were sufli(‘ient to 
provide very effective cooling (ref. 1). E. Schmidt (ref. 2), 
who proposed this cooling method, used coolant passages 


which were drilled into the turbine blade in a spanwise direc- 
tion in such a manner that the ends of the passages near the 
blade tip were closed whereas the inner ends near the blade 
root were open to a licpiid pool. When the blades are heated 
by the combustion gases and rotated, the heat flow into the 
cooling fluid which fills the blade i)assages creates tempera- 
ture differences within the coolant and in this way sets up 
rather intense free-convection currents. The Grashof num- 
ber which characterizes this flow is very high because of the 
large centrifugal forces, and the flow is consequently in the 
turbulent range. 

Very little experimental information exists on free- 
convection flow in the Grashof number range above 1 0*^ (refs. 3 
and 4). A theory developed in reference 5 to obtain informa- 
tion on the heat-transfer coefficients connected with this type 
of flow and on the thickness of the heated boundary layer, 
which builds up along the surfaces of the passage, starts from 
the assumption that the flow in a hole of small length-to- 
diameter ratio is the same as that for a vertical flat plate in an 
infinitely extended fluid. In the development of this theory, 
the shapes of the temperature and velocity profiles within 
the heated boundary layer were assumed. Sufficient experi- 
mental data could not be found in the literature to check the 
assumed profiles. 

The experimental investigation described in this report was 
conducted at the XACA Lewis laboratory in order to check 
these assumptions indirectly by obtaining a measurement of 
local heat-transfer coefficients in an arrangement which sim- 
ulates conditions that exist in the (coolant passage of a rotat- 
ing turbine blade. Difficulties were encountered in check- 
ing the velocity and temperature profiles because of unstable 
conditions in the test apparatus which will lx* dis(*ussed in 
RESULTS AND DISCUSSION. 

Detailed measurements in a rotating passage are practi- 
cally impossible because of the small dimensions involved. 
Therefore, a stationary apparatus, in which the centrifugal 
field which generates the free-convection flow in the rotating- 
blade coolant passage is replaced by the gravitational field, 
was designed, fabricated, and tested. Conditions for heat 
transfer should be comparable in both cases as long as the 
influence of Coriolis forces is neglected (ref. 5) and the 
Trandtl number for the fluid, as well as the Grashof numlx'r, 
is the same for the rotating and stationary conditions. 

This investigation was conducted at a Prandtl number of 
approximatelv 0.7 and in a Grashof number range from 10^ to 
10 ^^ 


1 Supersedes NACA RM E52F30, “Experimental Investigation of Free-Con vcction Heat Transfer in Vertical Tube at Large Grashof N urn I)crs,“ by E. R. O. Eckert and A. J. Diaguila, 
1952. 
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APPARATUS 

Basic considerations. -Exp('rimonts undortakon to study 
frce-convoction flow in a liquid-coolod turbine should be 
made in the same Grashof number range as, and with 
Prandtl numbers similar to, those of the cooling fluids used 
in the application (especially water and MlL-F-r)()24A 
grade JP-3 fuel). 

With regard to the Prandtl numbers, the following fact 
must !)(' eonsider(‘d: Gsually the results of experiments on 
free-eonveetion flow are |)lotted as Xusselt number against 
the product of Grashof and Prandtl numbers with the as- 
sumption that Xusselt number depends only on the product 
of both parameters. Kecent experiments and ealeulations 
(refs. 4 and 0 to 8) showed that this relation does not hold for 
a large l^randtl number range. Therefore', the Prandtl nuni- 
bei- of the test fluid should aj)|)roximate the l^randtl number 
of the fluid to be used in the application. 

The Grashof numbers which characterize the free-convection 
flow in a liepiid-cooled turbine are very large because of tlu' 
large centrifugal forces. They were found to extend from 
10*^ to lO’*^ (ref. 1). Special provisions are nec'essary to 
obtain such Grashof numbers in a stationary apparatus sinc(‘ 
the gravitational force is much smaller. The Grashof num- 
ber for free-convection flow under the influence of the gravi- 
tational force (j is 

0 >- ( 1 ) 

(Symbols are defined in the appendix.) 

A high Grashof number can be obtained in an experimental 
apparatus by use of a fluid with a large ratio, by providing 
an apparatus with sufficiently large chai-acteristic dimension 
/, or by a large temperature difference. The Pi*andtl num- 
bers 7V, the kinematic viscosities v, the thermal expansion 
coefficients p, and the /3/r J'atios ])resented in table 1 foj- a 
number of licpiids and gases are based on properties listed in 
reference 7 and other soui'ces. Mercury has the largest 13 /v~ 
ratio shown, but the use of this liquid in an apparatus with 
large dimensions is not convenient. In addition, mei-curv 
has a Prandtl number which is very low as compared with 
the Prandtl num.bei- of water or fuel, which are usually 
considered foi* turbine licjuid-cooling applications. Water 
has a ratio which increases with tempei*atui*e. However, 
the temperatures in an experimental setup would have to be 
ke])t well below 212° F to avoid evaporation when the 
apparatus is not pressui’ized. hlthylene glycol, oil, and sim- 
ilar licpiids have low /3/r^ ratios because of their high viscosity 
and are not well suited for the contem])lated experiments. 
Very large values of can be obtained w4ien licpiids are 
near their critical state. Values of for water and waiter 
vapor at a pressure of 200 atmospheres is plotted against 
temperature in figure' 1. Operation neai‘ the ci-itical state, 
however, has the disad vantage that the pi'operty value's 
change very raj^ielly. 

In general, gase^s have lowur values e)f p/v~ than liejuiels; how'- 
ever, the kinematic vise*osity, w Inch is inversely pro})ortional 
to the elensity, can be elee*reaseel anel 13/p^ ine'reaseel for gase's 



FiCiCre 1. — \'aliH's of ^/u- for water and water vapor lu'ar critical state 

(200 atm). 

by an increase in pre'ssure Air at a j)i*essure of 125 |)ounels 
per sejuare inch absolute, for instance, has a value coinpai'able 
te) that e)f water. In aelelition, the Praneltl number of air is 
in the same range as the Pi-andtl number of .\riL-P-5f)24A 
grade JP-8 fuel and w'ater at temperatures encountert'd in 
turbine litiiiid cooling. 
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A Top cover of configuration 1 
B Screening 
C Inlet 
D Outlet 

E Openings in steam jacket 
F Steam jacket 
G Condensate chamber 
H Exposed portion of condensate line 
J Top cover for configuration 2 
K Bottom plate 
L Burette board 
M Burette 

N Sj)herical condensate collector 
O Petcock 

P Heated-wall thermocouple 
Q Unheated-wall thermocouple 
R Steam thermocouj)le 
S Thermocou})le probes 

(a) Apparatus 

Figure 2. — Free-convection tube. 

From the previous considerations and for convenience, it 
was decided to use pressurized air in the apparatus. The 
apparatus consisted essentially of a vertical tube, closed at the 
bottom; the walls were heated so that heat was transferred 
to the pressurized air which filled the tube. Some provision 


had to be made at the top of the tube to replace the lieated 
air, which rose along the tube wall, with fresh cool air. The 
dimensions of the apparatus had to be (piite large in order to 
attain the desired Grashof numbei- j-aiige. Steam was used 
to heat the tube walls in order to ensure a uniform tempera- 
ture of the walls and to make the determination of tJie local 
heat flows convenient. The heat flow was based on the 
amount of condensate collected in small chambers. 

Tube arrangement. -The apparatus which was develo})ed 
from these basic considerations is shown in figure 2(a). 
Various component paiTs are designated by letter on the figure 
to simplify references in the text, ddie over-all dimensions 
of the tube ariu height, 13K feet; diameter, 24 inches; and 
wall thickness, % inch. The tube was fabricated from mild 
steel and coated with zinc on the inside and outside surfaces. 
The inside surface was ground smooth after the coating was 
applied. The height of the remaining roughness of the wall 
along the axial length adjacent to the steam chambers was 
measured and found not to exceed ±0.007 inch. This vari- 
ation was assumed to be ])ermissible on the basis of friction 
data obtained with forced flow th rough a tube with rough 
surfaces. 

Ail- at a temperature of approximately 80° F and a pressure 
of 23 to 125 ])ounds piu* scpiare inch absolute was introduced 
at the top cover A. Monel-wire screening B, which produced 
a large pressure drop, was installed at the air inlet to ensui’e 
uniform downward flow of the air. The wire screen con- 
sisted of two layers of 30X28 mesh (0.01 3-in. -diam. wire) and 
two layers of 50X40 mesh (0.0059-in. -diam. wire) stacked 
at 45° to each other. The desired pressure in the tube was 
obtained with throttle valves in the inlet and outlet lines (C 
and D, respectively). An orifice in the air-supply line was 
used to measure the air flow to the tube. 

Odie surface of the tube was heated with low-pressure steam 
su|)erheated by throttling to 2° or 3° F above the saturation 
tiunperature; the pressure ranged from 1 to 3 inches of mer- 
cury above atmospheric ])ressure. Four openings E along 
the length of the insulated steam jacket surrounding the tube 
were used to distribute the steam uniformly throughout the 
jacket F. A steam trap on the main condensate line and a 
throttle valve on the steam-su])ply line were used to obtain 
the desired steam pressure. 

Sixteen condensate chambers G, 7 % inches wide, were 
arranged along the length of the heated section of the tube; 
therefore, only part of the circumference of the tube was 
covered. The chambers varied in length from 6 inches at the 
top of the tube to 8)2 inches at the bottom and trapped only 
the condensate which developed on the section of the tube 
wall enclosed by these chambers. From each chamber, %- 
inch-diameter lines, shown in figure 2(b), carried the con- 
densate through the steam jacket to the condensate measuring 
apparatus. The exposed portion of these lines (H in fig. 2(a)) 
was heavily insulated to reduce heat losses. 

After the desired information was obtained with the 
arrangement shown in figure 2(a), the method in which the 
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(b) Tube with steam jacket removed showing condensate chambers 
and condensate flow lines. 

FioriiE 2. — Concluded. Free-con vection tube. 


air was iiitrodiiccHl to the tube was altc'red in ordcM* to study 
its (dlVct on heat transfer. The top cover A was l•el)laced hy 
the redesigned cover J. This cover admitted air over the 
entire tul)e cross section, whereas in the first cover a riii^ 
section was j)rovided around the wall to se])arate the rising 
heated air from the incoming, fresh, cool air. The new 
cover contained screening B consisting of a Js-inch-thick 
Fiherglas mat hetweeii two )i-inch-mesh (23-gage wir(‘) 
galvanized screens to ensure uniform air velocities. Tin' 
bottom plate K of the first configuration was replacc'd by a 
cover identical to cover J to permit the inti*oduction or 
exhaust of air at the bottom of the tube; this cover also made 
possible a study of superimposed free and forced convection. 

In a third configuration, the top cover or both the top 
and bottom covers A and K were removed so that room air 
could circulate into or through the tube. In one test a 
screen identical to the one in cover J was ])laced across the 
bottom opening to reduce the air flow. 

Condensate measuring apparatus. The condensate from 
the condensate chambers was measured on a burette board 
L. Each chamber G was connected to a burette M graduated 
in tenths of a milliliter. A spherical condensate collector 
N, l% inches in diameter, and a petcock 0 arrangement at 


the toj) of the burette were necessary for the following reason: 
Direct connections between the burettes and the lines from 
the condensate chaml)er H would permit steam to fill the 
|)ortions of the burettes not filled with condensate and thus 
provide an exposed area for steam condensation which would 
vary with the amount of condensate in the burettes. 

The burettes were insulated by two layers of air between 
two Plexiglas windows at the front and an insulating board 
at the rear. Nevertheless, the condensation caused by the 
heat losses from the i)urettes proved inconvenient in the 
first runs made with direct connections between burettes 
and condensate lines. The s])herical container N was 
emptied into the burette as soon as the condensate level had 
risen to about }i inch above the level indicated in figure 2(a) ; 
this ])rovided a relatively large volume of condensate and a 
small variation in exposed surface area on whi(*h steam would 
condense. The tube, inserted in the container N to slightly 
l)elow the center, provided a minimum condensate level in 
the sphere which could be conveniently reproduced. The 
time required to obtain this amount of condensate ranged 
from 1 5 to 30 minutes depending upon the operating condi- 
tion and the position of the condensate chambers. 

Instrumentation. -Five iron-constantan thermocouples P, 
si)aced as shown in figure 2(a), were used to measure the 
wall temperatures. These thermocouples were embedded 
in holes drilled from the outside into the tube wall to a 
depth of )4 inch. The thermocouj)le ends were welded into 
a small copper cylinder which fitted the holes closely. 
Thermocouple lead wires were attached to the outside wall 
of the heated tube and brought out of the steam jacket 
through the bottom flange. Fiberglas insulation and plastic 
coating protected the thermocouple wires from the steam. 
Two thermocouiiles Q were also used on the tube wall below 
the heated section to determine tube-wall temperatures. 
Stenm temperatures were measured by four iron-constantan 
thermocouples R sj)aced along the length of the steam jacket. 
Each thermocouple was so ai*ranged in a }^-inch-dianieter 

r- -Nickel wire, --SleeMubing, Copper wire, 

^ p.0004-in. diam. 0.25-in. O.D. 0.020-in. diam. 



'--Inconel wire. Ceramic \ 

003l-in. diam , , insulator--' 

( a) 


^--Butt-weld -Steel tubing -Ceramic 

thermocouple / 0.25-in. O.D. insulator 

junction 



(a) Resistance thennoiueter. 

(})) Iron-constantan thermocouple. 

Fku ue 3. — .Measuring probes for air-temperature survey. 
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tube that the junction was formed by welding the ends into 
a cap closing the tube. The thermocouple junctions were 
located midway between the tube wall and the steam-jacket 
wall. 

Originalh", the temperatures of the aii* inside the tube 
were to be measured with a resistance thermometer (fig. 3(a)). 
A bridge, a galvanometer, and a battery circuit, which 
measure the change of resistance with temperatures of a 
0.004-inch-diameter nickel wire, are used with this type 
of thermometer. Fluctuations of the air temperature, 
however, made reading with this instrument very difficult. 
It was therefore used only in certain runs to measure the 
air temperatures along the center line of the tube. 

In most cases the air temperatures inside the tube were 
measured with the thermocou])le probes shown in figure 3(b). 
The probe locations and identification numbers along the 
length of the tube ai*e shown by S in figure 2(a). The 
probes traversed from the center of the tube to the tube wall 
adjacent to the condensate chambers by means of a motor- 
driving arrangement on the probe mounting. With this 
arrangement it was possible to obtain temjierature readings 
of the heated air in five planes at any desired radius. Air 
temperatures were also measured with iron-(‘onstantan 
thermocouples in the inlet and outlet lines. 

The air pressure in the tube was measured with calibrated 
Bourdon type gages through pressure taps at two locations, 
one near the top and one near the bottom of the tube wall. 

Accuracy of measurements. — Calculation of the heat- 
transfer coefficients on the inside wall of the tube is based on 
the air temperature along the axis of the tube and on the 
temperature of the heated-wall surface. Since the tempera- 
ture near the center of the wall, where the thermocouples are 
located, is higher than at the inside wall surface, a calculation 
of the temperature drop due to the heat conduction through 
the wall was made; a temperature difference of approxi- 
mately 1° F was indicated. The accuracy of the potenti- 
ometer, the calibration accuracy, and the reading accuracy 
together added up to an error in the wall-temperature 
measurement of dz2.0° F. The calculated temperature 
difference between mid wall temperature and surface tem- 
peratures was within this range and was therefore disre- 
garded. The difference between the temperatures indicated 
by the thermocouples in the steam chest and in the wall was 
between 0.5° and 1° F. Steam temperatures and air 
temperatures were estimated to be accurate within ±1.5° F. 
Repeated readings of the tempei*atm*e of the heated air at a 
specified location, however, were only within ±3° to ±5° F 
depending on the location of the thermocouples. These 
fluctuations of the air temperature will l)e discussed in 
RESULTS AND DISCUSSIONS. 

The mercury U-tube measurements of the steam pressure 
were accurate within ±0.05 pound per square inch. The 
saturation temperature determined from this pressure agrees 
within 1° F with the steam temperature measured by the ther- 
mocouples. Air pressures were measured with two Bourdon 
type gages. The difference between the readings of the 
gages was small, and the average of both was used to deter- 
mine the air density. The accuracy of these measurements 


determined by calibration was found to be 2 percent below 
30 pounds per square inch absolute and 1 percent above this 
value. The condensate was collected during the runs and 
measured volumetrically. Duration of the runs was deter- 
mined with a stop watch. The measurements of the rate of 
condensation for the complete run could be reproduced 
within ±2.0 percent. The flow orifice for the supply air was 
calibrated and estimated to be accurate within 1 percent. 


CALCULATION OF LOCAL HEAT-TRANSFER COEFFICIENTS 


EQUATION DEFINING HEAT-TRANSFER COEFFICIENT 


The local film heat-transfer coefficients inside the tube at 
various locations along its length were determined from the 
following equation: 


Th 


. Qc 


( 2 ) 


In this equation the area A considered for heat transfer was 
that section of the inside surface of the tube which was en- 
closed by the individual condensate chamber. Included in 
this value is one-half the thickness of the walls of the con- 
densate chambers attached to the tube wall. The total 
thickness of the condensate-chamber wall is 0.05 inch. The 
temperature of the air ta in equation (2) was obtained along 
the tube axis directly opposite the center of the condensate 
chamber. Only five air-temperature measurements were 
made along the axis; therefore, the air temperatures opposite 
the centers of the individual condensate chambers were 
determined from a cuiwe faired through these measured 
temperatures. The temperature measured in the center 
of the tube wall was used in equation (2) as the wall tem- 
perature t,c. This temperature was constant along the tube. 


EVALUATION OF CONVECTIVE HEAT FLOW AND HEAT LOSS 

Determination of convective heat flow Qc. — The total heat 
How Q was determined by measurement of the condensate 
collected in the individual burettes. This condensate, how- 
ever, consisted not only of the amount collected in the indi- 
vidual condensate chambei*s (G in fig. 2(a)) but also of the 
amounts formed by the heat flow to the surrounding atmos- 
phere from the part of the condensate lines outside the steam 
jacket (H in fig. 2(a)) and from the spherical condensate 
collectors (N in fig. 2(a)). In addition, heat is transferred 
from the interior tube surface not only to the air by convec- 
tion but also to the cooler top and bottom sections by radia- 
tion. Correspondingly the convective heat flow per unit 
time Qc Ri^iy be determined from the e(j nation 


Qc— Q— Qa — Qr (3) 

where 

Q = wh (4) 

The heat flow to the surrounding atmosphere from the con- 
densate lines and in the spherical collectoi*s is denoted by 
Qa and the radiative heat flow by Qr in equation (3). In 
equation (4), the condensate weight flow per unit time is w, 
the total measured rate for the individual burette correspond- 
ing to the condensate chamber (or location along the tube 
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surfac(') undrr coiisidc'ration. The Icm’IU h is tlu' heal of 
vaporization |)or |)oimd of waU'i* at saturation iirossuiv, that 
is, tho inoasuird prossuro in the steam ja(*kot. The small 
dogreo of siipi'rhoating was nogloctod. The heat flow Q 
obtained by use of ecjuation (4) is therefore the total heat 
released by condensation of tlie weight w. 

Determination of heat loss Qa - Qr. In order to deter- 
mine the convective heat flow Qc, it is necessary to subtract 
from the total heat flow the heat lost to the atmosphere Qa 
from the exposed condensate lines and collectors, and the 
heat i-adiation Qr. Both quantities were determined in 
heat-loss runs. If the tube could be evacuated, heat would 
be removed from the interior tube wall only by radiation 
towards the cooler top and bottom sections. The conden- 
sate collected in such a run was the result of this radiative 
heat flow and the heat losses in the condensate lines and 
sj)herical collectors Qa. 

Actually, the laboratory facilities which were connected to 
this apparatus did not provide complete evacuation of the 
tube. Tlu'refore, a small convective heat transfer still took 
place in the interior of the tube, and equation (3) holds for 
such heat-loss runs as well as for heat-transfer runs. The 
convective hcait flow is, however, much smaller because of 
the lower density of the air and because the heated air inside 
the tube is not replaced in these runs but is heated to a tein- 
perature near that of the wall. This convective heat flow 
(Micountered in heat-loss runs may be indicated by Q[ and the 
amount of condensate collected in such a run by w' . Then 
equation (3) becomes 

Q: = w'}i-Q,-Qr (5) 

When the i-oom conditions are the same, the heat lost from 
the condensate lines and the spherical collectors Qa and from 
radiative flow Qr is the same as that in the actual heat- 
transfiM* lains; tliis heat loss can be comjiuted by use of e((ua- 
tion (o). 

If e([uations (3) to (5) are combined, 

Q^=wh — v/li + Qc (h) 

In this e([uation w is the amount of condensate collected 
during a heat-transfer run, w' is the amount collected during 
a heat-loss run, and Ql is the convective heat transfer in a 
heat-loss run. The value of Ql is small compai-ed with the 
other terms in ecpiation (h) under the test conditions set uj) 
in the heat-loss runs and can be determined by an ajiprox- 
imate calculation. The equation for the average Nusselt 
number in free-(‘onvection laminar flow (ref. 4) on a vertical 
plate is 

A^=0.o55(6V/V)i/^ (7) 

The local Nusselt number which follows was obtained from 
this relation by multiplying by the factor 3/4 (I’ef. 7, p. 162): 

A%=0.416(6V/V)'/4 (8) 

The laminar equation was used to obtain the convective 
heat-transfer coeflicient for the conditions existing in the 


tube during the heat-loss runs, because the low density of 
the air in the tube and the small temperature difference 
lowered the Grashof number for these runs below the critical 
value. The convective heat flow Qi was calculated by use 
of equation (2) and the temjieratures t,^ and t^ which were 
measured during these runs. All data were then availal)le 
to calculate the convective heat flow Qc in equation (6). 
The convective heat flow Ql encountered during the heat-loss 
runs amounted to 2 to 4 percent of the convective heat flow 
Qc in the heat-transfer runs. 

Alternate method of determining heat loss. -A calculation 
was also made to obtain an estimate of the heat flow by 
radiation Qr from the heated walls to the unheated toj) and 
bottom sections of the tube. In this calculation, the tem- 
perature of the tube surface was assumed to be 212® F and 
that of the bottom circular surface, 90® F. This temjierature 
range is rejiresentative of the severest conditions. The 
emissivity of the galvanized iron surface was taken as 0.20. 
For the heat-transfer runs, the heat exchange by radiation 
between the center of the wall section enclosed by the lowest 
condensate chamber and the l)ottom ])late was (*alculated on 
these assunqitions to be apiiroximately 0.75 ])ercent of the 
convective heat flow Qc. At other locations along the tube 
wall, the radiative heat transfer was still smaller; therefore, 
the heat transfer by I’adiation could be neglected. 

The heat loss Qa in the condensate lines and the spherical 
collectors can be determined by arranging tra])s at the loca- 
tion where the condensate lines leave the steam jacket. 
These tra])s were so designed that no condensate coming 
from the condensate chambers could flow into the s])herical 
collectors. These collectors then (‘ollected only condensate 
formed on their surfaces and that which formed within 
the condensate lines outside the steam jacket. From the 
condensate measured in the burettes by this method, the 
heat loss Qa could be found. The trajis were arranged 
only on a selected number of condensate lines in order to 
check the heat losses obtained by the first method described. 
(Convective heat flows Qc determined in this manner agreed 
with the values obtained by the first method within 3 percent. 

EXPERIMENTAL PROCEDURE 

Heat-loss runs. — The ])ur])Ose of the heat-loss runs was to 
determine the heat loss Qa (in the condensate lines and the 
s])herical collectors) and the radiated heat Qr, as described 
in the section Determination of heat loss Qa-I Qr, and also 
to obtain the information necessary for calculating the con- 
vective heat flow Q'^. In the heat-loss runs, the tube was 
])ai‘tially evacuated by reducing the air pressure to 7 inches 
of mercury absolute, the minimum permitted by the available 
laboratory altitude-exhaust facilities. I^efore the data were 
lecorded, the tube was heated by the low-pressure steam for 
aiiproximately 2 hours. This length of time was required 
for the tube-wall tenqieratures of the heated and unheated 
sections and for the temperature of the air in the tube to 
reach equilibrium. After eciuilibrium was obtained, the 
(luautity of condensate in the burettes was recorded at 20- 
minute intervals, and an average of at least five readings 
was obtained. The heat -loss runs were repeated at various 
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Figure 4. — Typical air-tcnipcratiire distribution alonp; tube axis and 
wall in heat -loss run. 

inUM'vals Ihroiighout llie investigation. Air temperatures 
in the tube ami tul)e-wall temperatures in a heat-loss run 
are presented in figure 4. From these temperatures and the 
Nusselt numbers given 1)V equation (8), it was ])ossil)le to 
calculate the heat transfer Q' in the tube during the heat- 
loss runs. The heat losses Qa tind Qr determined in this 
manner increased from 20 percent of the total heat-flow 
rate Q at the high-i)ressure runs to 45 ])ercent for the low- 
pressure runs. 

Heat-transfer runs. -Heat-transfer runs were made for a 
series of air pressures inside the tube between atmospheric 
and maximum air-sui)i)ly ])ressure (125 Ib/sq in. abs) in 
order to obtain the desired range of Grashof numbers. 
Runs were also made with different amounts of air flowing 
through the upper part of the tube. The ])urpose of this 
variation was to detect the influence of air-llow rate on the 
air-temi)erature distribution in the tube and on the heat 
transfer. For each run, the tube was preheated for a period 
of approximately 2 hours and data were then taken at 20- 
minute intervals over a period of 2 to 8 hours. Throughout 
each run, conditions were carefully kept constant. The 
meiisurements produced all necessary data for calculating 
local heat-transfer coefficients. In addition, the temperature 
field of the air inside the tube was investigated by measure- 
ment of the temperature piofiles along the cylinder radii 
in the planes indicated in figure 2(a) by the probes S. 

In the course of the investigation, the method by which 
the heated air in the tube was replaced was changed in ordei- 
to study its effect on the temperature field and the heat 
transfer. For this purpose, the top and bottom covers 
were changed as described in the APPARATUS section. 

RESULTS AND DISCUSSION 

Temperature field. -The air-temperature distribution 
within the tul)e is shown in figures 5 to 9. The results of 


measurements on the first cover configuration, shown in 
figure 2(a), are presented in figures 5 and 0. The data 
obtained with the second cover configuration are presented 
in figures 7 and 8. The air-temperature measurements of 
other configurations studied are shown in figure 9. 

Temperature profiles along the tube radii measured at 
different tube heights are i)resented in figures 5 and 7; 
figures 6, 8, and 9 show temperatures measured along the 
tube axis. Included in these figures are inlet-air and wall 
temperatures. Temperature profiles measured with probe 1 
are not shown because the temperature fluctuations in this 
plane were too large to permit reproducible measurements. 
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Figure 5. — Radial air-temperature profiles at four probe locations 
along axis for first cover configuration. .\ir pressure, 85 pounds per 
scpiare inch absolute; air flow, 0.2.50 pound per second. 
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(a) Hums with air ])rossurc‘s from 125 to 04 j)ouncls per scjuare inch 

absolute. 




(b) Huns with air pressures from 44 to 25 pounds |)er square inch 

absolute. 

Ficjuke 0. — Temperature of inlet air, air along tube axis, and wall 
for first cover configuration. 

AVheii this investigation was started, it was assumed that 
the air which is heated along the cylindrical surface of the 
tube would rise within a thin layer along this surface and 
would he removed through outlets D (fig. 2(a)) and that the 
air temperature in the central core of the tube would he 
(fuite uniform. Figure 5 shows that the essential tempera- 
ture drop actually occurs very near to the tube wall. How- 
ever, figure 6 indicates that an appreciable increase in tem- 
perature with distance from the air inlet occurs in the central 


core; the amount of increase is approximately 30° F. In- 
spection of the various curves reveals that the differences 
l)etween the individual curves in figures 6(a) and (b) cannot 
be connected with the pressure of the air in the tube nor 
with the amount of air passing through the top portion of 
the tube. Such a tem])erature distribution (with warmer, 
and therefore lighter, air layers below cooler, and therefore 
heavier, ones) creates a very unstable state which can exist 
only in the ])resence of dynamic and viscous forces connected 
with flow. It is therefore understandable that the ffow in 
the central core is fluctuating and turl)ulent. This condition 
is also indicated by the temperature fluctuations encountered. 
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Figure 7. — Radial air-temiierature profiles at four probe locations 
along axis for second cover configuration. Air pressure, 99 pounds 
per square inch ab.solut(‘; air flow, 0.508 pound p(‘r second. 
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PTgtre 8. — Temperature of inlet air, air along tube axis, and Avail for 
second cover configuration. 


especially in the upper poidion of the tube. These lluet na- 
tions are probably responsible for the heat transferred from 
the hot boundary layers alon^ the wall into the central 
portion of the tube. 

A (piestion then arises as to what degree the temp(‘iatui‘e 
distribution in the core of the air is influenced by the specific 
method used in removing the warm air stream along the 
wall from the top of the tube and replacing it by fresh cool air. 
One arrangement investigated involved re])lacing the top 
cover A (fig. 2(a)) by cover J. The bottom plate K was 
exchanged for a cover identical to the top cover J; however, 
no air Avas admitted or removed through the bottom cover, 
and thus its action was not different from that of plate K. 
Air was admitted through the top cover J and removed 
through theexhaust lines D. Theresultsof the air-temperature 
measurements are presented in figures 7 and 8. Again, 
a thin heated boundary layer is indicated near the tube wall. 
Also, a temperature increase with distance from the air inlet 
along the tube axis is present, and the amount of the increase 
throughout the tube is approximately 30° F. The various 
curves, which are closer together than those of figure 0, 
indicate more stable conditions than for the first configura- 
tion. The temperature measurements in the central plane 
(fig. 8) were made with the resistance thermometer, whereas 
all the other temperature measurements were made with 
iron-const ant an thermocouples. 

A more radical change in the air-flow conditions was ob- 
tained in the following manner: The top cover was removed 
so that the tube was opened to room air. The hot boundary 
layer discharged into the room was thus replaced by cool 
air flowing into the top of the tube. The temperatures 
along the tube axis, which are shown in figure 9 (curve A), 
indicate an increase with distance from the inlet of approxi- 



Figure 9. — Temperature of inlet air, air along tube axis, and Avail for 
various cover configurations. 

mately the same magnitude as before. When part of the air 
is admitted into the bottom of the tube through screening 
consisting of a }^-inch-thick Fiberglas mat between tAvo 
)i-inch-mesh (23-gage wire) galvanized screens, a tempera- 
ture increase Avith distance from the lower inlet occurs in 
the loAver portion of the tube, whereas the temperature field 
in the upper portion is essentially the same as it was with 
the previous configuration. Curve B (fig. 9) shows the 
temperatures for this case. Only Avhen the screen at the 
bottom was removed and the tube AA^as open at both ends 
did the temperature along the axis become nearly uniform 
(fig. 9, cuiwe C). A temperature distribution similar to 
curve B Avas obtained in the second cover configuration when 
part of the air Avas admitted through the upper and part 
through the Ioaatm* cover (fig. 9, curve D). 

The investigation indicates that a temperature variation 
along the tube axis is consistently associated Avith the specific 
configuration Avhich provides free-convection floAv in a tube 
closed at the bottom and heated along its Avails. Such a 
configuration corresponds to the liquid-cooling arrangement 
for turbines proposed by Schmidt. The magnitude of the 
temperature increase is approximately the same regardless 
of the cover configuration. In addition, the essential differ- 
ence in the temperature distributions for the individual runs 
is a general change in the temperature level which may be 
determined by the degree of mixing of the outgoing and in- 
coming air streams. 

In figure 10, the ratio of the temperature difference between 
the wall and the air along the center line of the tube to the 
same temperature difference near the inlet (probe 1, fig. 2(a)) 
is plotted against the tube length for the first cover con- 
figuration, the second cover configuration, and Avith the top 
cover removed. Xo apparent difference in this temperature 
ratio exists for these configurations. For tlie Grashof and 
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Figure 10. — Ratio of temperature differences along tube length be- 
tween wall and air at center line to same temperature difference nt‘ar 
inlet. 


J^raiidtl iuiinl)(*r ranges investigated, tlie curve laired through 
the observed points represents the average temperature 
increase of air in a vertical tube closed at the bottom, heated 
along the wall, and having a leiigth-to-diameter ratio of 5. 

In order to obtain additional information on the air flow in 
the core of the tube, a scale model one-sixth the actual size 
was fabricated from pyrex glass with a top cover similar to 
J shown in figure 2(a). The model was not heated because 
the object was to obtain only information on how the flow 
in the tube is affected by the inlet and outlet configuration. 
The model was filled with smoke and then clear air was intro- 
duced at the top. A photograj)h of this apparatus during the 
smoke test is shown in figure 11. It may be observed that 
the smoke is quite dense in the lower third of the tube, indi- 
cating that very little mixing of the air occurs in this portion 
of the tube. In the central third of the tube, the smoke is 
still present, but it is less dense, indicating a slight dqgri'c of 
mixing. In the upper third of the tube, the air is completely 
clear as a result of considerable mixing with the air supplied 
at the top. The Reynolds number for the air flow under the 
conditions shown in figure 11 was lower than in the full-scale 
test apparatus. Visual observations, however, were made 
up to Reynolds numbers which almost equaled those in the 
full-scale appaj’atus. \o appi’cciable change in flow pattern 
could be observed for increasing Re.ynolds number. 

Heat-transfer coefficients, — Local heat-transfer coeffi- 
cients were calculated by use of equation (2). Local coeffi- 
cients were converted to dimensionless Xusselt numbers; a 
characteristic length I, which will be discussed later in this 
section, was used in this conversion. Xusselt numbers are 
plotted in figure 12 against the product of Grashof and 
Prandtl numbers. This parameter is customarily used in 
free-convection correlations and it was retained here al- 
though, as pointed out in the APPARATUS section, recent 



Figure 11. — Smoke test of flow inside glass scale model of tube 
indicating penetration of incoming air by decrea.se in smoke intensity. 

investigations indicate that it does not apply for a large 
Prandtl numlier range. 

The property values of the air were introduced into the 
dimensionless parameters at a reference temperature equal 
to the arithmetic mean between wall temperature and air 
temperature at the tube axis. The readings on the burette 
connected with the highest condensate chain bei* were con- 
sistently low. Consequently, the coefficients for the section 
of the tube adjacent to the highest condensate chamber were 
omitted in figure 12. The range of Grashof numbers extends 
from 10^ to 10^^, which is higher than the ranges covered in 
previous investigations. In these figures, curves representing 
local Xusselt number relations for free-convection flow on a 
vertical flat plate in the laminar and turbulent ranges are 
also included. The line representing heat transfer for 
laminar flow corresponds to equation (8). The following 
relation for average XTisselt number is given in reference 5 
for free-convection heat transfer under turbulent flow 
conditions: 

l^=0.0210(GrPry'^ 

This relation, which is in agreement with relations given in 
references 3 and 4 for the range of this investigation, shows 
that for turbulent flow the Jieat-transfer coefficient is 



Nusselt number, Nu Nusselt number, Nu 
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(a) First cover configuration. 



Product of Groshof and Prondtl numbers, GrPr 


(b) Second cover configuration. 

Fiot’RE 12. — Dinnuisionless correlation of local free-convection heat-transfer coefficients in closed tube and comparison with known equa- 
tions for vertical flat plate. 
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Product of Groshof and Prondtl numbers, GrPr 
(c) \’arioiis cover configurations. 

Figi'RE 12. — Conclu(l(‘cl. Dimensionless correlation of local free-convection Iieat-transfer coefficients in closed tube and comparison with known 

equations for vertical flat plate. 


])roportioiial to the 1 /5 power of the (listanee from the lower- 
plate end and that the local Niisselt niimher is 1.20 times 
the average value. Therefore, the relation for the local 
heat-transfer coefficient is 


2/.5 


A%=0.0252(G'/7T) 


These equations hold for a vertical plate with uniform tem- 
perature, but the wall temperature of the tube was constant 
only where it was surrounded by the steam chest (fig. 2(a)). 
In the portion below the steam chest, the wall temperature 
decreased with distance from the air inlet. In order to 
account for this fact, a fictive end of the heated section was 
established as indicated in figure 13. In this figure, the 
measured wall temperatures are plotted against distance 
from the lower end of the steam chest. The air temperatuj*e 
along the tube axis as extrapolated from the measured values 
is also indicated on the figure. It can he seen that the tem- 
perature difference between the wall and the air decreases 
with distance below the end of the steam chest and reaches 
a value of zero 8 inches below the jacket for this particular 


run. It can be expected that in a first approximation con- 
ditions are similar on a tube with a constant temperature 
difference which has its lower end at a position such that the 
two dashed ai’eas in figuj*e 13 ai’e eciual. In this way, the 
expression f which is important for the heat 

transferred, becomes ecjual for the real tube and the one witli 
a constant temperature difference. The fictive end of the 
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Fig r RE 13. — Determination of fictive lower end of steam chest for 
comparison of local heat-transfer data with flat-plate values. 
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steam chest was determined in this way. The chaj*acteristic 
length I contained in the Nusselt and Grashof numhers was 
measured from the fictive end of the steam chest to the center 
of the condensate chamber being considered. This length 
influences Nusselt and Grashof numbers only for the two 
lower condensate chambers. Figure 12(a) shows the Nusselt 
and Grashof numbers calculated for two points, with the 
effective length L measured from the lower end of the steam 
chest to the centers of the lower two condensate chambers. 
These points are connected by arrows, with corresponding 
points calculated by use of the distance from the fictive end 
of the steam chest foi* the characteristic length. 

For all arrangements investigated, figure 12 indicates that 
the Nusselt numbers obtained generally correlate satisfac- 
torily as functions of the Grashof number. In the laminar 
range, the measured values are above the curve for the local 
laminar Nusselt relation for a vertical flat plate. In the 
turbulent range, the mean value of the data is displaced below 
the curve for the local Nusselt relation for a vertical Hat plate 
by approximate!}" 35 percent. The fact that the degree of 
scatter of the measured values for the second cover con- 
figuration (fig. 12(b)) is less than for the first configuration 
(fig. 12(a)) confirms the conclusion that the flow conditions 
in the core were more stable for this configuration. 

The figures also show that the way in which the heated 
ah* is removed and replaced by fresh air does not influence 
the heat-transfer coefficients when these are determined by 
use of local differences between wall temperature and air 
temperature along the tube axis. The experiments indicate 
that it is possible to obtain an approximate value of the heat- 
transfer coefficients in tubes closed at the lower end and 
heated around the surface for a small length-to-diameter 
ratio (approx. 5) from known relations for vertical flat plates. 
A calculation of the heat transferred requires the knowledge' 
of the air temperature along the tube axis. An approximate 
value of the temperature increase of the air along the axis 
can be obtained from the curve in figure 10. 


SUMMARY OF RESULTS 

The following results were obtained in an investigation of 
free-convection flow within a stationary tube closed at the 
bottom, heated around the periphery, and having a h'ligth- 
to-diameter ratio of approximately 5; the heated air was 
continuously replaced by fresh cool air at the top. Condi- 
tions that exist within the coolant passages of a i-otating 
turbine blade were simulated in the tube. 

1. The greatest drop in the temperature of the air in the 
tube occurred very neai- to the tube wall. This drop in- 
dicated that a thin boundary layer existed along the heated- 
wall surface. 

2. The temperature of the air in the core of the tube in- 
creased with distance from the air inlet. Such a temperature 
distribution creates unstable flow conditions. The magni- 
tude of the temperature increase did not appear to be in- 
fluenced by the specific air inlet and outlet configurations. 

3. The method of removing and replacing the heated air 
did not inffuence the heat-transfer coefficients when these 
were determined by means of local temperature diff’erences 
between the heated vv-all and the air along the tube axis. 

4. Heat-transfer coefficients, which were converted to 
dimensionless Nusselt numbers, correlated as functions of the 
Grashof numbers. The range of Grashof numbers investi- 
gated extended from 10® to 10^®, which is beyond the range 
previously investigated . 

5. In the turbulent range investigated, the mean valuesof the 
Nusselt numbers were approximated 35 percent below known 
relations for vertical flat plates. These relations may be 
used to obtain approximate heat-transfer coefficients in 
heated tubes having a length-to-diameter ratio of approx- 
imately 5. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 30, 1952 
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APPENDIX 

SYMBOI.S 


'Hie following symbols are used in ibis rej^ort: 
inside area of section of tube surrounded by condensate 
chamber, s(j ft 

specific heat of fluid at constant pressure, Btu/(lb)(°F) 
(jrrashof number, « 

V 

acceleration due to gravity, 82.2 ft/sec^ 

local heat-transfer coefficient, Btu/(°b") (s(i ft) (sec) 

average heat-tj*ansfer coefficient, Btu/(°F)(s([ ft) (sec) 

heat of vaporization of water, Btu/lb 

thermal conductivity of fluid, lUu/(°F) (ft) (sec) 

characteristic length, ft 

local Xusselt numbei*, IHIk 

average Xusselt number, 1 1! Ik 

Prandtl number, 

heat flow, fku/sec 


/ (empi'ratui’e, °F 

w condensate weight (low, Ib/sec 

P coefficient of thermal expansion of fluid, 1 /°R 
jjL dynamic viscosity of fluid, lb/(sec) (ft ) 
p density of fluid, Ib/cii ft 
V kinematic viscosity of (luid, p/p, s(i ft/sec 

Subscripts; 

.1 atmos|)her(‘ 

a ail’ 

c convectiv(‘ 

R radiation 

tv wa 1 1 

Superscript: 

' values obtained in heat-loss runs 
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TABLE I.— PROPERTY VALUES CONSIDERED IN FREE- 
CONVECTION STUDIES OF VARIOUS FLUIDS 


1 

Prandtl 

number, 

Pr 

1 

Kinematic 
viscosity, v, 
S(i ft/.sec 

1 

1/°K 1 

1 ! 

sec*2 

(°RT('ft)^ 

Liquids at tenipriatun' ol i»8° K (saturated stale) » 

Mercury 

0. 0249 

0. 12:iX10-s 

1 

^ 0. 101X10-3 

007X 10' 

M 1 F-5(i24 A ( J P-3) fuel 

12.3 

1.04 

. 725 

07. 1 

Water 

7.02 

l.OH 

. 10 

8.5 

Water (212° F) 

1.74 

.310 

.35 

:351 

Ethyleneglycol 

204 

20. 0 

. 30 

. 0845 

Engine oil (unused) 

10. 400 

970 

.39 

.0000415 

Gases at temperature of 122' 

F and pressure of 14.22 pounds per square inch “ 

Freon-12 

0. K4 

3. 17X10-' 

1.72X10-3 

17. 1X10' 

Carbon dioxide 

.79 

10. HO 

1.72 

1 . 40 

Ammonia 

.91 

18. 9 

1.72 

.483 

Air 

. 70 

19. 89 

1.72 

.434 

Air (125 Ib/sq in.) 

.70 

2. 20 

1.72 

33. 65 

Hydrogen (212° F') 

. 05 

174. 2 

1.49 

.00492 


“ Excej)t where noted. 

This value in error in ref. 7. 
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